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Abstract. The FUSE far-ultraviolet interstellar spectra toward seven targets in NGC 6231 show that the molecules H2, HD, 
and CO as well as various atomic species are distributed in several clouds. The main absorption component found on the sight 
lines lies in the Lupus cloud region at a distance of about 150 pc, and there is a weaker second one, presumably in the Sco OBI 
shell surrounding NGC 6231 (d 1.8 kpc). H2 excitation modelling is used to constrain the radiation field and the density in 
the gas; HD is used to estimate the abundance of H + . The small angular separation of the target stars allows column density 
variations to be probed over the field of view, on scales of 0.05 pc in the case of the Lupus cloud distance. They are 40% for 
H2 and 60% for H 1. A rather strong radiation field inside the molecular clouds suggests a separation of the gas into smaller 
cloudlets also along the line of sight. 

Key words. ISM: structure - ISM: molecules - ISM: abundances - structure - Ultraviolet: ISM 



1. Introduction 

When observing the interstellar medium (ISM) using absorp- 
tion line spectroscopy, one has to live with the restriction that 
the data gives just pencil beam information along single lines 
of sight towards selected background sources. Absorption line 
analysis will never provide as wealthy structural information as 
the analysis of emission lines or imaging. 

However, one can get closer to this goal by choosing groups 
of close together background sources, in this way enlarging the 
information density per angular element. 

Open star clusters, rich in young UV bright background 
sources, play an important role in such projects when dealing 
with far-UV (FUV) absorption line spectroscopy of the ISM. 
In particular, the open cluster NGC 623 1 in the Sco OB 1 as- 
sociation provides at least seven early-type stars in its central 
region which are suitable as background sources for FUV spec- 
troscopy, and it still shows moderate extinction in this wave- 
length region for a far-away object in the Milky Way disk. 
The cluster is located in the Galactic plane (Z = 343.46°, b = 
+ 1.17°), at a distance of ~ 1.8 kpc (Raboud et al. 1997). The 
E(B - V) values for its stars vary between 0.43 and 0.46 mag 
(Baume et al. 1999, see Table 1). 

We will analyse the sight lines towards these seven back- 
ground sources with respect to absorption by the Galactic in- 
terstellar medium, paying particular attention to the variation 
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in the physical parameters of the gas over small angular scales. 
The background sources span a total area of 5.8' x 3.1'. At the 
distance of the cluster this means a linear extent of the gas of 
about 3.0 pc xl .6 pc. Part of the absorbing gas will very proba- 
bly be located more nearby, resulting in an even higher spatial 
resolution. 

The spectral range of the FUSE instrument covers the 
Lyman and Werner bands of molecular hydrogen H2 and of 
deuterated hydrogen HD, the H 1 Lyman series from Ly /? on 
bluewards, the CO E-X band, and various metals. We are thus 
able to derive physical parameters for the gas like excitation 
temperatures and densities from excitation and column density 
ratios. Sect. 2 introduces the data themselves and the analysis 
methods used; Sect. 3 gives an overview of the results and a 
discussion of their reliability. Sect. 4 contains the interpreta- 
tion, followed by concluding remarks in Sect. 5. 

2. Instrumentation, observation, and data 
reduction 

2.1. Data observation and preparation 

The FUSE instrument consists of four coaligned telescopes 
with Rowland spectrographs and two microchannel plate de- 
tectors. The spectral range of operation is between 905 and 
1187 A, at a resolution of about A/AA =* 15000 or 20 kms" 1 
(FWHM). An extended overview of the instrumentation and 
the on-orbit performance of the satellite is found in Moos et al. 
(2000) and in Sahnow et al. (2000). 

The open cluster NGC 623 1 was selected as a target with 
the criteria in mind of a) a cluster with a number of hot early 
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Table 1. List of the targets and observational parameters. The spectrum of HD 152270 was observed in 5 separate pointings, thus no unique ID 
is given. The E(B - V) values are from Baume et al. (1999). 



Target 


Sp.Type 


RA (2000) 


Dec (2000) 


FUSE ID 


Exp. date 


Exp. time [s] 


E(B - V) 


HD 152200 


09.5III 


16 53 51.69 


-41 50 33.1 


B0250101000 


2001-08-09 17:48:00 


5367 


0.43 


HD 152219 


09.5III 


16 53 55.65 


-41 52 51.8 


B0250201000 


2001-08-10 13:46:00 


4905 


0.46 


HD 152234 


B0.5Ia 


16 54 01.88 


-41 48 23.4 


B0250301000 


2001-08-11 18:22:00 


5123 


0.45 


HD 152249 


09.5Iab 


16 54 11.66 


-41 50 57.8 


B0250401000 


2001-08-12 02:26:00 


4947 


0.44 


HD 326329 


09V 


16 54 14.06 


-41 50 09.0 


B0250501000 


2001-08-09 12:45:00 


5485 


0.44 


CPD -47 7712 


BOIV 


16 54 00.44 


-41 52 44.0 


B0250601000 


2001-08-09 04:24:00 


9611 


0.46 


HD 152270 


WC7 


16 54 19.70 


-41 49 11.5 


B0 13020x000 


2001-07-08 23:00:48 


2950 


0.46 




1 6 n 55 m 00 



20 s 54 m 0C' 5 40 s 
a (2000) 



20 s 53 m 00 s 



Fig. 1. Positions of the target stars. The background shows the 60 fim 
IRAS emission map (see intensity scale at right). 



type stars (earlier than about B3, see Table 1) providing clean 
UV continua without perturbing stellar atmospheric absorption 
and b) a moderate colour excess so that the UV flux is not di- 
minished by a too strong interstellar grain absorption. Table 1 
lists the FUSE targets and the observational parameters. As 
an additional benefit the individual background sources are all 
separated by at least 1 ', thus we are able to use the large LWRS 
aperture (30" x 30") for observation. The separation is large 
enough to assure that always only one UV bright star is located 
within the aperture. A map of the NGC 6231 region and the 
locations of the FUSE targets is given in Fig. 1. 

We use the CalFUSE standard reduction package (version 
2.0.5) for the basic data reduction and calibration. The different 
exposures for each target are cross correlated and coadded for 
each detector channel. The coadded spectra are then rebinned 
to a binsize of 4, which results in a resolution of 0.028 A per 
bin (or 8.4 kms -1 at 1000 A), still well below the instrument's 
spectral resolution. 



2.2. Fitting procedure 

The analysis of the absorption lines is only possible if one 
applies full Voigt profile fitting to the absorption as a whole. 
Strong damping wings and heavy blends especially of the 



molecular lines hamper the identification of single absorption 
lines and the determination of equivalent widths for them; 
blends of lines of different species or excitation levels are com- 
mon. 

The Doppler parameter of the absorption b was estimated 
beforehand using a standard curve of growth analysis, but on a 
restricted set of absorption lines. This was done separately for 
the molecular and for the atomic gas since they possibly exist 
in different environments. For the molecules we used the H2 
transitions of the J3, J4, and J5 excitation states which, for our 
sight lines, define the full flat part of the curve of growth (as 
preparatory plots showed). These lines thus tend to have only 
moderate damping wings, if at all, and readily provide equiva- 
lent widths. For the metals we used the full series of available 
iron absorption lines except the few that are obviously blended. 

The b values so estimated were then taken as input param- 
eters for an interactive^ 2 fit to the FUSE spectra. We generate 
a full theoretical Voigt profile absorption line spectrum over 
the FUSE wavelength range, fitting the radial velocity of the 
absorbing gas and the column densities of all major absorbing 
species (including Hi) at once. A continuum level is estimated 
by eye. This way we can also handle the extended damping 
wings and blends in the lines in particular of the lower H2 ex- 
citation levels and also the various inter-species blends. In the 
process of fitting the initially estimated b values were adjusted 
as well, if nessessary. Fig. 3 shows a sample of curves of growth 
using the final b values and column densities. 

The lines' transition parameters are taken from Abgrall 
et al. (1993a,b) for the H2 Lyman and Werner band 
(utilising the electronic version of the data provided by 
S.R.McCandliss 1 ), Morton (1975, 1978) and Wright & 
Morton (1979) for HD, and Morton & Noreau (1994) and 
Federman et al. (2001) for CO. The metal transitions are from 
Morton (2003), with some values for Fe taken from Howk et al. 
(2000). The full set of transitions is included in the fitting pro- 
cess, however the weakest lines do not neccessarily contribute 
to the result. A small number of stronger transitions had to be 
disregarded due to unidentified blends or detector effects. 

We use the wavelength/velocity calibration that is provided 
by the FUSE lALiF spectra. Although the FUSE wavelength 
calibration is, at present, known to be problematic (see, e.g., 
the discussion in Par. 3.4 of Wakker et al. 2003) comparison 
with the velocity structure observed in high resolution optical 
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Fig. 2. Absorption lines from the seven target star spectra. The sample portions are binned (to ~ 8.4 km s -1 ), but not normalised to a continuum. 
The An line at 1048 A is located in the left wing of a strong H 2 absorption feature. Stellar wind features affect the continuum in the WR 
star HD 152270. For some sight lines (HD 152249, HD 152270, HD 326329) a second blueward velocity component is clearly visible, while 

it is much weaker in the others. The dashes above the spectra mark the positions of the main velocity component at 8 kms~' and of the 

weaker second component at -40 km s -1 . At the left end of the box of the Ni line the transition of Ni (1 134.4 A) can be seen, shortwards of 
the HDR(0) line a blend of the H 2 transitions R(3), 4-0 and R(6), 5-0. 



spectra (Crawford 1995, 2001) shows a high similarity. We thus 
assume that the wavelength scale defined in the lALiF (and 
also the 2ALiF) spectra is very close to reality. Corrections 
would have had to be applied to the spectra of some of the 
other detectors, though, but they are not used in this work. 

Fig. 2 shows samples of various unblended absorption lines 
on the seven sight lines. In Fig. 4 a sample of typical heavily 
blended H2 absorption structures is given, together with the full 
Voigt profile fit. 

2.3. Error estimates 

Errors for the column densities are determined as 3 cr devia- 
tions from the fit. In addition to these, we have to take care 
of additional uncertainties from four sources which are not in- 
cluded in the fitting process itself. 

The FWHM: 

The instrumental FWHM in our spectra seems to change 
over the FUSE detector segments, with the lALiF spectra 
having a slightly higher FWHM (24 kms~', in the case of 
HD152200) than the 2ALiF spectra (22 kms" 1 ). This adds 
some uncertainty to the column density fit, however it mea- 
surably affects only the fits for N 1 (near A — 1 1 34 A) and for 
the J = 3 to 5 levels of H2. In general, the variation of column 
density with FWHM is in the range of the small error from the 
X 2 fit. For the few elements mentioned one may get up to three 
times this value. 

The three H2 levels are the ones which we use in fixing the 
b value because of their location on the flat part of the curve 
of growth. The curve of growth procedure, however, deals with 
equivalent widths only and is thus independent from changes in 



FWHM. Thus, the origin of this particular variation on chang- 
ing FWHM is rather caused by the general uncertainty in fitting 
elements on the flat part of the curve of growth, not by the un- 
certainty in FWHM itself. 

The lines of Ni are in general difficult to fit due to their 
close triplet nature, they are highly sensitive to changes in any 
fitting parameter, including FWHM. 

The continuum: 

Another question is the dependency of the resulting col- 
umn densities on the assumed continuum. There seems to be 
no proper way to quantify this dependency. Of course, elements 
with only few transitions are more sensitive to changes in the 
local continuum of a particular line. Among those, the saturated 
lines of especially 1 and N 1 are more affected, since the fit is 
here mainly sensitive to changes in the flanks of the absorption 
profile, while the profiles themselves are not very sensitive to 
changes in column density. 

Overall, we find a variation in column density of about the 
size of the^- 2 fit errors (about twice that value for 1 and N 1) if 
one varies the local continuum by about 10% of the originally 
assumed level. The resulting continua are outside the range of 
what would be considered a valid continuum for our spectra; 
therefore this gives a good estimate of the uncertainty intro- 
duced by the continuum selection. 

The b value: 

The b value is the largest and also the easiest to underesti- 
mate contributor of errors, whether one uses the conventional 
curve of growth technique or any sophisticated fitting proce- 
dure. Even in high-quality spectra it can be difficult to fix. 

For the metals, only iron is suitable for determining b, all 
other elements having too few usable transitions in the FUSE 
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Fig. 3. Sample curves of growth. Curves of growth were only used to determine an initial preliminary b value for the fitting process; this b value 
was then subject to some adjustment during the fit. Shown are single-component curves of growth for the b values that resulted from the fit, but 
the equivalent widths plotted are sums over the full absorption features, i.e. sums over the equivalents widths of the individual blended velocity 
components. For this reason a discrepancy between equivalent widths and the curve of growth is visible for the curve of the two-component 
sight line HD 152270. 



spectral range. As a consequence, some of the latter are also 
more error-prone. In our case, the column densities in particular 
of oxygen and nitrogen are extremely sensitive to changes in b, 
with their line strengths at the level of the flat part of the curve 
of growth. The less abundant elements are located on the linear 
part of the curve and are thus well defined, but for the same 
reason they cannot be used to fix the b value. 

The various H2 excitation states are located all over the 
curve of growth, but only the slightly saturated middle levels 
define the b value for the molecular gas. As for the metals, their 
column densities again are most sensitive to changes in the b 
value, while the equivalent widths of the transitions and their 
profiles change only marginally. The quality of the fit therefore 
does not suffer much and does not provide any measure for the 
quality of the b value. 

Fitting the spectrum with different b values gives an im- 
pression of the variation. We varied the originally determined b 
value by + 1 km s _1 . Column densities for elements having only 
absorption lines with a strength of the flat part of the curve of 
growth changed by up to 0.5 dex, for elements low on the linear 
part or high in the damping part about 0. 1 dex or even less. As 
error from the b value uncertainty we assume, for each element, 
the mean deviation for^ 2 fits with b + 1 km s _1 . 

The stellar contribution: 

Stellar atmospheric features are relatively few in number 
and then easily visible, mostly wide features. Strong wind pro- 
files are found, especially in the Wolf-Rayet star HD 152270. 
We tested modelling the stellar spectra using Kurucz atmo- 
spheric models, without a satisfying result. Broad stellar fea- 
tures are therefore included by eye in the assumed continuum. 



If interstellar lines are affected by narrow stellar atmospheric 
lines in a way that such a continuum cannot unambiguously be 
defined, these regions are excluded from the fit. 

The individual errors from all these sources -\ 2 fit and un- 
certainties in FWHM, continuum, and b value - add up to the 
given total errors. In most cases, the errors from the fit itself 
are very small, so the larger contribution comes from the un- 
certainty in the continuum, in the FWHM, and especially the b 
value. 

3. Discussion 

The column densities resulting from the fit are listed in Table 2. 
If a species is not detected we give 3cr upper limits determined 
at the positions of strong, but not neccessarily the strongest, 
unblended transitions. 

3. 1 . Velocity structure 

We find one prominent absorption component, and in some 
spectra a second one bluewards which is weaker but clearly 
distinguishable. 

The first, stronger velocity component is seen in H 2 absorp- 
tion in all spectra at heliocentric velocities of about -8 kms~'. 
Velocities can be fitted with a precision of 1 km s _1 or below. 

We observed that the velocity of some metal absorption 
lines is different from that of the molecules by up to 6 kms~' 
bluewards. This is observable only in those spectra which 
show a single velocity component. In the other spectra it is 
probably obscured by blendings of the components. We fitted 



O. Marggraf et al.: Intermediate scale structure of the ISM towards NGC 6231 



5 




1050 1052 
Wavelength [A] 



Fig. 4. Sample portion of the HD 152234 
spectrum with a one-component fit, in 
which the positions of the fitted transitions 
are indicated. The continuum count rate 
in this region of the spectrum is ca. 6000 
(binned). The full wealth of H 2 transitions 
mainly of the (4-0) vibrational states is vis- 
ible. Also visible is a strong HD line of the 
J = state (HD R(0), 4-0). For the J = 1 
state of HD we measure only upper lim- 
its. The An transition at 1048 A (see also 
Fig. 2) is affected by some unknown absorp- 
tion on its left wing, maybe by a weak sec- 
ond component bluewards. Also other lines 
hint at the presence of a weak second com- 
ponent. 



these lines using slightly different velocity offsets. In two cases 
the Oi transition at 1039 A is seen to be shifted by another 
~ -4kms _1 . 

The origin of this shift may be unresolved substructure in 
the form of an additional blueward component that is visible 
only in metal lines. A complex velocity substructure is known 
to exist from ultra-high resolution optical spectra of sight lines 
in this general direction (Crawford 1995, 2001). A problem in 
the FUSE wavelength calibration would not be restricted to the 
transitions of metals and thus can be ruled out as a cause for 
such a shift. 

The second, less prominent velocity component is found 
in at least three of the spectra (HD 152249, HD 152270, 
HD 326329), mainly in the weaker metal transitions at a ve- 
locity of about -40 kms~'. It should be noted that these three 
sight lines are located closely together at the eastern edge of 
the target field. On other sight lines a second component may 
also be present at that velocity judged from well-defined metal 
lines like Pn(1152). 

In optical ultra-high resolution spectra, Crawford (2001) 
resolves numerous velocity components in Nai and Ki, with 
stronger column density peaks around -39 kms -1 , -9 kms -1 , 
and +2 kms" 1 . Crawford (1995) finds numerous absorptions 
of CH, CH + , and CN, again peaking around velocities of 
-9 kms -1 and +2 kms -1 . Especially the velocity of the 
-9 kms 1 component varies, however, depending on the en- 
vironment the absorbing element lives in. 

CO emission line surveys for this general direction show 
two main velocitiy components as well (e.g., Bronfman et al. 
1989). However, the component velocity separation in that 
data is slightly larger than in our spectra (about 40 kms -1 
compared to our 32 kms~'), and the velocities are shifted by 
-10 kms -1 . As supposed by Bronfman et al. (1989), the emis- 
sion at negative velocities should originate in the far back- 
ground of NGC 623 1 . It is therefore not related to the absorp- 
tion in the NGC 623 1 star spectra, despite the similar velocity. 



The emission components seen at positive velocities can be as- 
sociated with the Lupus molecular cloud region at a distance of 
about 150 pc (Crawford 2000), and those may indeed interfere 
with our sight lines. 

3.2. Column densities 
3.2.1. Molecules 

H2 can be fitted with high accuracy due to the multitude of tran- 
sitions covering a large range of transition probabilities. Only 
the middle excitation states (esp. J = 2,3) are less accurate, 
resulting from their location on the flat part of the curve of 
growth (see above, Sect. 2). A mean NQI2) - 20.38 indicates 
a sizeable amount, in line with the values of E(B - V). The 
higher excitation levels J — 6 and 7 which are highly desired 
for excitation modelling are hardly separated from the noise, 
though. 

Also the HD absorption is, in the 7 = state, strong enough 
for a clean fit. The J = 1 state is close to the detection limit but 
may be visible in some cases. We give upper limits for the HD 
J — 1 column density. 

Fitting CO is more difficult. The resolution of FUSE is too 
low to resolve the individual transitions in the CO absorption 
bands. Typically, the Doppler b value in an environment where 
CO is found is low. This means that the equivalent widths 
are not neccessarily located on the linear part of the curve of 
growth; they are thus no longer additive. The column densities 
can therefore only be lower limits. 

The CO abundances along our sight lines are, with respect 
to the H2 abundances, within the usual range. CO column den- 
sities in the ISM typically increase significantly at H2 column 
densities > 10 19 cm 2 . The ratio of N(CO)/N(H 2 ) found fits 
into the relation drawn by Federman et al. (1980). However, 
the ratio of JV(CO)/N(H + 2H 2 ) deviates from other Galactic 
sight lines, presumably due to the fact that the total sight line 
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Table 2. List of column densities log /V [cm -2 ] towards NGC 6231. FWHM, radial velocities v, and Doppler widths b are in kms~'. Upper 
limits are at least 3cr. The b values are "effective" values for the measured absorption lines and should be lower for individual unresolved 
subcomponents. 





152200 


152219 


152234 


152249 




152270 




326329 




-417712 


FWHM 


24 


24 


24 


22 


22 


22 


22 


22 


22 


22 


V 


-8 


-8 


-8 


-8 


-40 


-8 


-40 


-8 


-35 


-10 


b 


7.0 


7.0 


7.0 


6.0 


6.0 


8.0 




7.0 


7.0 


7.0 


H 2 , J = 
H 2 , J = 1 
H 2 , J = 2 
H 2 , J = 3 
H 2 , J = 4 
H 2 , J = 5 
H 2 , J = 6 
H 2 , J = 7 


20 14+ 002 

1 Q Q7+0.02 

18 12 + °- 2 Z 

1 7 40+ - 45 
i ; ,4U_ 45 

15 96 +016 

15.40!°;!° 

< 14.15 


20. io:™ 2 

20.00!°;° 2 
18 39 + !H! 

lo,J7 -0.11 

17 56 +0 - 43 
i ; . ju „ 43 

15 91 +018 
1J - 71 -0.18 
1 c i a +o m 
±J.lt_ 007 

< 14.00 

14 06 +0 02 


20 2r 02 

_ Q03 

20 01 +0 02 

_ 002 

18 25 +016 

17.18^| 
15 94+016 

1 c- /-»^4-0 OS 

15.22!™ 8 
14 09+° 05 
14 12 +002 


20 00 +0 02 
z-u.uo _ 005 

19 89 +0 02 

_ 002 
18 42 + °°, 8 

10 '^ -0.08 

15.88t!J 

1 ^ nr+0(M 

1 4 95 

< 14.00 

< 14.20 


> 16.00 
14 70 + ° 32 
1 4 Q0+ - 26 
14 30+ 017 
14.20 + °;| 3 

- 


20 22 +0 - 02 
20 03 +0 - 02 
18 68 +0 05 

io.uo _ 005 

17.84!^ 
15 65 +012 

1 i /Ai-i+O 03 

14 -99^o4 
< 14.15 
14 02 +0 04 


- 
- 

- 


20.08^ 
19.95 + °;° 2 

18 58 + !!!5 

1 7 Q7+0.36 
1 '- yl -0.36 

< 15 96 +0 3 ° 

i c i o+0 07 

O.i8_ ; 08 

< 14.15 

< 14.10 


< 15 40 +0 49 

< 14.60 + ° 32 
15 12 +006 

1J ' iz -0.06 

< 14 48+0 25 

1 * * ^4-0 48 

< l4 - 4 6-048 

- 


20 25 +0 - 02 

20 ir 002 

18 72 +0 - 06 

° -0.06 

18.26^11 
16 29 +014 

1 5 7/ 
1J,i -0.09 

< 14.20 

< 14.20 


HD, y = 

hd, y = i 


14.64^ 
< 14.00 


1 4 59+0-06 

< 13.90 


14 62+ 05 

-0.04 

< 13.95 


14 5r 006 

< 13.90 


- 


14 55+ - 07 
< 14.15 




14 6r 010 

< 14.20 


- 


14 7 j +0.09 
1H - ' x -0.10 
< 14.30 


CO 1 
CO 2 


13.85:°* 
> 13.70 


i t no +0.08 
U./8_ 004 

> 13.63 


13 91 +0 03 
i j. ?! _ 002 

> 13.69 


13.82^ 
> 13.67 


- 


13 95 +0 - 05 

1J - 7J -0.02 

> 13.60 




13 99+° 11 
> 13.79 


- 


13-89!^ 
> 13.87 


V 


-6 


-2 


-6 


-1 


-40 


+2 


-40 


-3 


-40 


-4 




9.0 


9.0 


10.0 


10.0 


- 


9.0 


- 


10.0 


- 


9.0 


Ci( 3 P ) 
Ci( 3 P 2 ) 


15 15 +0 - 05 
14.86!- 3 
< 14.60 


15 05+ 03 
u.uj _ 006 

< 14.90 

< 14.60 


15 17 +002 

1J - i '-0.02 

15 00 +0 03 
< 14.75 


15 05 +0 07 

< 14.75 

< 14.45 




15 00+ 015 

iJ ' uu -0.15 

< 14.80 

< 14.60 




15 14 +005 
1 4 97+O.O6 
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'from fit, 2 from W A (CO C-X @ 1088 A); see Sect. 3.2.1 



samples large amounts of Hi in regions where no molecular 
gas, and especially no CO, exists. 

The molecule column densities, of H2 as well as of HD 
and CO, are not automatically linearly related to those of H 1, 
because they have their origin in a far more constrained fraction 
of the sight line. Molecular gas is in general confined in the 
denser regions of the interstellar clouds and is thus of a more 
clumpy nature. 

3.2.2. Neutral hydrogen 

The deviation of the column densities of H 1 is of prime im- 
portance. They are the basis for the depletion values and thus 



any analysis of their variation, as well as for any statement on a 
possible abundance variation between the different sight lines. 
Judged from the relative abundances of the, in general, only 
marginally depleted P and of Fe (see, e.g., the review by de 
Boer et al. 1987), all our Hi abundances are consistent with 
the typical P and Fe values observed along high column den- 
sity sight lines (Fig. 5). The values towards HD 152234 and 
HD 152249 fully agree with the results from the IUE survey of 
Diplas & Savage (1994) as well. 

We observe a variation in the total gas column density, 
N tot = N(H 1) +2/V(H 2 ), between log Mot = 21.34 (HD 152200) 
and 21.52 (HD 152270), which is about 50% over an angle of 
5f7 (60% if one only takes Hi into account). The variation is 
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Fig. 5. Element depletions as derived from 
the seven sight lines observed. The width 
of the error boxes varies for better readabil- 
ity only. The bold lines represent the values 
typically found in the interstellar medium in 
warm (upper bold line) and cool (lower bold 
line) gas (Savage & Sembach 1996). The 
scatter especially in Ar, O, and Si is due to 
difficulties in determining the column den- 
sities for these metals. 



rather a linear increase from the western to the eastern edge of 
the field of view than a fluctuation due to structure on even 
smaller scales. The expected large number of Hi absorbers 
along the sight lines will, even if they are spatially structured, 
smear out any fluctuation, leaving just a general smooth in- 
crease of column density. This increase in N tot is in agreement 
with the map from IRAS for this region, which shows a smooth 
rise in dust emission in the same general direction (Fig. 1). The 
effect may also be connected to the second velocity compo- 
nent that we observe towards the group of targets located in the 
north-eastern region of our field. 

3.2.3. Metals 

The column densities derived for the metals are given in 
Table 2. Fe n is the only metal with a large number of usable 
transitions in the FUSE spectral range, which results in a reli- 
able fit and the smallest uncertainties in column density. 

Also Pn seems, with only two transitions (AA 1124, 
1152 A), to give good results if one compares its derived de- 
pletion with typical values (de Boer et al. 1987). To get proper 
metal depletions we use the metal column densities summed 
over all velocity components, since we are unable to separate 
any velocity components in Hi. This, however, reduces the 
value of the depletions for the interpretation of the physical 
conditions in a particular velocity component. 

The C i ground state transitions blend with various transi- 
tions of the excited states of neutral carbon and are difficult 
to unravel. Therefore, upper limits are listed in most cases. 
Information on C n, which is the dominant state of carbon in 
neutral gas, is not available in the FUSE spectral range. 

For some important ions the number of usable transitions in 
the FUSE spectral range is small (Ni AA 1134.165, 1134.415, 
1134.980, 1159.817; O i A 1039.230; Sin A 1020.699), and we 
cannot rule out stellar contamination. We therefore have to as- 
sume large errors for these column densities. 

The accuracy of the An abundance is compromised due to 
difficulties in fitting the only two transitions. The second veloc- 
ity component at -40 kms" 1 seems to be present on all sight 
lines for An, which has an effect on the b value to choose for 
the fit. Also, the An lines are located on the flat portion of 



the curve of growth, making a column density determination 
highly sensitive to the chosen b value. 

Oi is especially hard to fit on the two-component sight 
lines, in particular towards HD 152249 and HD 152270. Fits 
are possible only to one flank of the profile, which requires a 
precise wavelength calibration at the position of the transition. 
Table 2 states the results from the fit for O i, but the errors given 
are only the results from the fitting process itself and cannot 
represent the larger true uncertainty. 

4. Interpretation 

4.1. Sight line structure 

NGC 623 1 is located in the centre of the Sco OB 1 association, 
at the near side of the Sagittarius spiral arm. The structure of 
the interstellar medium along the line of sight is complex. 

The extinction survey of Neckel & Klare (1980, field 
343/+ 1) shows in this direction a sharp increase in interstellar 
extinction by A v = 0.75 mag that takes place between about 
100 and 300 pc distance. As mentioned before (Sect. 3.1), 
the sight line passes the outer regions of the Lupus molecu- 
lar cloud region, an extended star forming region at the near 
edge of the Sco-Cen OB association. The 12 CO emission maps 
of Bronfman et al. (1989) show that our sight line misses this 
CO cloud by 1°, only about 2.6 pc at a distance of 150 pc. 

Another sharp increase in extinction is found by Neckel 
& Klare (1980) at a distance of 700-1000 pc, again by A v = 
0.75 mag. Beyond this point, the sight line seems to be clear up 
the cluster at about 1900 pc. 

The heliocentric radial velocity of NGC 6231 is about 
-25 kms -1 , depending on the method used in the determina- 
tion (e.g., Laval 1972; Dambis et al. 2001). Gas surrounding 
the association is expected to have been cleared away, so its 
absorption velocity will be more negative. In consequence, we 
can assume that gas at more positive absorption velocities is 
not associated with the cluster but located somewhere on the 
line of sight. 

Following this reasoning, the velocity component seen at 
-40 kms"' probably originates in the outwards-moving shell 
of Sco OBI. Ultra-high resolution optical spectra towards 
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Table 3. H 2 Boltzmann temperatures r , and equivalent excitation 
temperatures T 35 for the -8 km s _1 component. 



Target 


T 0i 


[K] 


^35 


[K] 


HD 152200 


66 


+ 2 


309 


+ 


13 


HD 152219 


71 


+ 3 


259 


+ 


11 


HD 152234 


64 


+ 2 


286 


+ 


11 


HD 152249 


70 


±5 


235 


+ 


6 


HD 152270 


65 


+ 4 


267 


+ 


41 


HD 326329 


69 


+ 3 


216 


+ 


1 


CPD-417712 


68 


+ 3 


204 


+ 


10 



HD 152249 (Crawford 2001) show substructure in this velocity 
component that is not resolved by FUSE. 

It seems reasonable to assume that the absorption at 
-8 kms -1 is associated with the outer regions of the Lupus 
molecular clouds where strong absorption is expected to hap- 
pen. However, velocities at this longitude do not provide, using 
Galactic rotation models, any distance information that would 
back up this assumption. The detection of the second veloc- 
ity component at v ra d = -40 kms -1 in only a subgroup of the 
sight lines may also point to a location of that component in the 
background, near Sco OB 1 . 

The slightly velocity-shifted metal lines then may result 
from an unresolved blend of absorption at the Lupus cloud 
distance and at larger distances, maybe including even more 
unresolved components not containing H2. The extinction ob- 
served between 700 and 1000 pc may be connected to these 
unresolved components. 

4.2. Physical parameters of the gas 

The excitation of the H2 molecule results from two separate 
mechanisms: The lower levels are predominantly excited by 
collisional excitation that is described by a Boltzmann temper- 
ature. The higher levels (from J = 2 on upwards) are popu- 
lated mainly by molecules cascading back from even higher 
excitation states after either being pumped to these states by 
UV photons, or in the molecule formation process. This ex- 
citation is in general described by fitting another Boltzmann 
distribution characterised by a so called equivalent excitation 
temperature. The low level H2 excitation temperatures for the 
individual sight lines are given in Table 3. 

We did extensive excitation modelling of H2 for some of 
the sight lines. Appendix A gives a more detailed description 
of the model. In short, the model assumes a single isothermal 
cloud of constant density and the equilibrium population of the 
rotational levels J = to J = 7 for a given total H2 column 
density is calculated. 

Using simplifications such as a constant interstellar radia- 
tion field and an approximate treatment of the UV attenuation 
by dust inside the cloud, the models allow an estimate of the 
level of the FUV radiation field and of the hydrogen density. 

For most stars it is not possible to fully fit the excitation 
using a single H2 component model. Interestingly, this may 
again point to a hidden two component structure in the gas, 
with physically distinct environments. 
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Fig. 6. Measured H 2 excitation and excitation model for HD 152234. 
Testing of various models restricts the parameters to a hydrogen den- 
sity n of not more than 1000 cirr 3 and a radiation field at A 1000 A 
with U A = (3.1 + 0.5) ■ 10~ 15 erg cm 3 A" 1 . Shown is the fit for 
n = 100 cnT 3 . 



As it was shown, e.g., by Flower et al. (1986), the shock 
front of an expanding shell can also excite H2 in a way that 
is comparable to relatively high Boltzmann excitation tempera- 
tures. There are, however, indications against a current shock in 
the Sco OBI shell (Crawford 2001) and also no indications for 
a shock in the Lupus clouds (Crawford 1995). The Nai/Can 
ratios indicate a previous shock in the Sco OBI environment 
that released part of the Ca bound in dust grains, but since then 
the gas has cooled down again. 

Gas in the Sco OB 1 shell, however, would have a different, 
more negative, velocity and thus would not hide in the main 
velocity component. The most plausible explanation for a hid- 
den second component is found in either substructure inside 
the Lupus cloud region or in the gas causing the strong extinc- 
tion at 700 to 1000 pc distance. From our data it is not possible 
to make any definitive statement on this or on any even more 
complicated velocity structure. 

The excitation on the HD 152234 sight line allows a fit of a 
one component model (Fig. 6). Fitting a range of models con- 
strains the gas density n to not more than 1000 cirr 3 and the 
radiation field to 1/^(1 000 A) = (3.1 ±0.5) • 10~ 15 erg cm- 3 A" 1 . 
This value is about 100 times the typical Galactic disk radia- 
tion field (e.g., Habing 1968), which could hint at a significant 
source of radiation closeby. 
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HD 152200 


HD 152219 


HD 152234 


HD 152249 


HD 152270 


HD 326329 


CPD-417712 


Mean 


E(B - V) 


0.43 


0.46 


0.45 


0.44 


0.46 


0.44 


0.46 




Z(H 2 ) 

S(HD) 

log/ 


20 37+ - 01 
14 i2 +am 
-0.66 
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c -7 1 +0.06 
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<; 74+0.03 

J - ' -0.02 

6 52+ 003 

-0.02 


5 67+ 005 

J '"'-0.05 

6 44+ 005 

U '+ H -0.01 


c -7 o +0.07 
3 ' /8 -0.06 

6 49+° 05 

u -^-0.03 


5 62+°-° 9 

J ' u -0.08 

6 34+o.n 

u -' -0.05 


5 70+°- 08 

J - ' -0.08 

6 61+°" 


5.70 + 0.05 
6.50 + 0.09 



Another estimate of the conditions in the gas is provided 
by the ratio of the two lowermost excitation levels of HD. In 
contrast to H2, transitions between even and odd rotational ex- 
citation states are possible. The population of the J = and 
7=1 states of HD therefore depends not only on collisional 
but also on excitation by the prevailing radiation field. 

However, if one assumes the temperature, density, and ra- 
diation field derived from the H2 excitation to be valid also 
for HD, one would get a much higher population of the HD 
J = 1 level than is observed (N(HDj =1 )/N(HD J=Q ) £ 0.45 
instead of the observed 0.2). The ratio that is observed con- 
strains the radiation field in the gas to not more than Ua - 
1.4 • 10~ 15 ergcirT 3 A -1 (or / < 35), which is at least a factor 2 
smaller than the value from the H2 models. 

To explain this discrepancy, one has to recall that the H2 ex- 
citation model assumes a single cloud. A high radiation field is 
needed to penetrate this large cloud and produce the observed 
excitation pattern. A number of smaller clouds with the same 
overall column density along the sight line would show an ex- 
citation pattern similar to the observed one already at a smaller 
radiation field, provided the radiation is not attenuated by other 
cloudlets. 

We tested additional models for a series of clouds with only 
a fraction of the total column density each. The observed ex- 
citation pattern is equally well fitted by two or four smaller 
clouds and by a radiation field smaller by the same factor, plus 
a smaller b value for each cloud. Such a model would still be 
consistent with the larger number of absorption lines and the 
correspondingly low b values Crawford (2001) finds from op- 
tical ultra-high resolution spectra. 

This suggests a considerably smaller radiation field than 
required by the simple one-cloud model (possibly not much 
higher than the Galactic radiation field) and a number of 
smaller distinct clouds along the sight line. 

The modelling of the H2 excitation furthermore allows us 
to estimate the formation rate of H2 in the gas. The H2 dissoc- 
itation rate from the one-component model is 1.2 • 10~ 13 s _1 . 
This would require a formation rate coefficient R ^ 1.6 • 
10~ 16 cm 3 s _1 (for a density of n = 100 ctrT 2 ) to maintain 
equilibrium between formation and dissociation of the H 2 . 
However, the commonly accepted value for the H2 formation 
rate is 1-3 • 10~ 17 cm 3 s _1 , depending mainly on the compo- 
sition of the interstellar dust. Provided that the H2 formation 
takes place at only the usual rate, the local radiation field will 
then dissociate the H2 gas within a time scale of roughly 10 6 a. 



A stable condition might exist if the gas density is as high as 
~ 500 ctrT 3 , thus lowering the required formation rate coeffi- 
cient. 

Crawford (1989) gives an analysis of the CH column densi- 
ties towards Sco OBI. Using a simple model with an assumed 
density of n > 100 ctrT 3 , he reproduces these column densities, 
but the model fails to reproduce the CH+ abundance. Provided 
that H2 and CH coexist in the same region of space, this con- 
strains the density in the molecular gas to be in the range of 
100cm 3 SbS 1000cm 3 . 

The quality of most of the metal abundances is too low for 
any constraints on physical gas parameters. 

4.3. Molecular abundances 

The ratios of the molecular species, H2/HD and H2/CO, are, 
within the range of the errors (Table 4), constant over the sight 
lines. This allows a good estimate of these astronomically im- 
portant parameters. 

We get a mean log /V(H 2 )/log /V(CO) of 6.50 + 0.09. 
Although the physical conditions are different, this value can 
be compared with results from other models, e.g. from van 
Dishoek & Black (1988, especially their model HI) for the 
f Oph sightline. As expected, this locates the NGC 6231 gas 
at the low density edge of these translucent cloud models, at 
a radiation level probably near the Galactic disk average (e.g., 
Draine 1978) and densities of about 500 crrr 3 . 

Also, we obtain a mean log A^(H 2 )/log /V(HD) of 5.70+0.05. 
This value is not related to the H/D ratio in a linear way but it 
also depends, due to the formation processes involved, on addi- 
tional parameters like the proton density in the gas. Also, in all 
dominant processes involved in the production of HD ionised 
hydrogen or deuterium plays a role. For this reason the ra- 
tio of H2/HD leads to an estimate for the ionisation level in 
the molecular cloud. Assuming a density of n — 100 cm -3 
and an optical depth at 1000 A of t uv - 2 we get a frac- 
tion of ionised hydrogen of x(H+) ^ 3 • 10~ 5 . We assume 
N(D)/N(H) = 1.2 • 10~ 5 , the mean value from recent measure- 
ments (Moos et al. 2002; Hoopes et al. 2003; Bluhm et al. 1999) 
on sight lines beyond the Local Bubble. This low ion fraction, 
however, does not allow conclusions on the abundance of free 
electrons in the gas, which are mainly provided by C n under 
the prevailing conditions. 
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4.4. Excited neutral carbon 

The excitation of neutral carbon and its dependence on gas 
density and temperature is described by de Boer & Morton 
(1974). The ratio of Ci( 3 Pi)/Ci( 3 P ) is probably < -0.2 dex 
for our gas, on the HD 152234 sight line roughly -0.2 dex. 
With this value and the temperature of 64 K that was derived 
above (Sect. 4.2) for the H2 gas, one gets a gas density of 
«h - 130 cm 3 . Alternatively, a gas density of 100 ctrT 3 as 
derived from the H2 models would require a gas temperature 
of T - 76 K for the observed C 1 excitation. In general, the gas 
pressure n H T from Ci is around 8000 K cm 3 . 

4.5. Depletions 

The depletion of the metals seems, with the uncertainties in 
mind, not unusual. The values for phosphor and iron are in 
the range that is typical for cool disk gas. For An one might, 
from the depletion diagram (Fig. 5), crudely estimate an addi- 
tional underabundance by -0.2 dex due to partial ionisation. 
This would correspond with an ionisation rate of the gas of 
0.07, which would require a photoionisation rate of T(H) = 
2.2 • lO" 16 s _1 (assuming T = 5000 K and n H = 0.1 ctrT 3 ) or, 
density independent, r(H)/n e = 3.1 • 10~ 14 s _1 . However, the 
uncertainty in the Ar 1 column density is large (as discussed in 
Sect. 3.2.3) and these rates should not be overinterpreted. 

4.6. Effects of velocity substructure 

In all these considerations we have to bear in mind that there 
is velocity substructure in our absorption lines which we are 
unable to resolve, and that there is probably also absorption at 
similar velocities from physically distinct regions in the ISM. 
The gas seen in the partially ionised An is surely distinct from 
the gas seen in the molecules, the Ar 1 being in the outer bound- 
aries of denser regions or possibly warmer intercloud medium. 
Also the number of molecular clouds is unclear, or its sepa- 
ration into smaller cloudlets along the line of sight. The CO 
is expected to exist in the shielded cores of these clouds, so 
its density will be somewhat higher than that of the H2 gas. 
Nevertheless, all methods and all models seem to point at an 
extended substructure of small cloudlets which can be deeply 
penetrated by the surrounding radiation field. 

4.7. Spatial substructure 

Besides the structure of the gas along the lines of sight, a main 
aspect of this paper is the structure over the field of view, per- 
pendicular to the sight lines. Of the metals only Fen and Pn 
provide useful information on this structure. To get a valid com- 
parison with the H 1 abundance that is not resolved in velocity, 
one has to use the column density sums over the velocity com- 
ponents observed, which results, again, in a lack of significant 
variation for both species. Thus, the statistical blend of different 
absorbers hides any information on column density fluctuations 
over a particular cloud. 

The molecular gas may, however, provide some insight 
since it very likely is more sparsely distributed along the sight 



lines. We observe a rather irregular variation of the H2 column 
density over the field of view, with variations of 40% over a 
scale of roughly 1' (HD 152219 and CPD-41 7712 being the 
extremes). 

The angular separation of these two sight lines corresponds 
to only ~ 0.05 pc at the distance of 150 pc assumed for the 
molecular gas. Under the - simplifying - assumption of circu- 
lar H2 cloudlets, one derives a radius of not more than 0.1 pc 
(23) for the strongly absorbing H2 cloud in front of CPD - 
41 7712, at the edge of the triangle formed by HD 152200, 
HD 152219, and HD 152249. 

The difference in column density in this region is about 9 • 
10 19 ctrT 2 . With a radius of 0.1 pc such a cloutlet would then 
have an average density of roughly 300 cirT 3 , again in accord 
with the values derived above. 

Some variation is also visible in the equivalent excitation 
temperatures that are connected to the intensity of the radiation 
field at the location of the gas. These temperatures, however, 
are only slightly anticorrelated with the H2 column density. 

The imperfections in the calibration of the FUSE spectra do 
not allow us to make any hard statement on potential velocity 
substructure over the field of view. From the agreement of the 
individual H2 absorption velocities we assume a quite homoge- 
neous velocity distribution of at least the H2 gas over our field 
of view. 

Absorbers seen in metal lines are far more abundant along 
the sight lines, so the absorption is seen at many more different 
velocities. The variation seen in the velocities of these lines 
probably rather comes from faint variations in the blends of 
different absorbers than from real velocity substructure. 

5. Concluding remarks 

We find a main absorption component at a velocity of 
-8 kms~', and a weak second one at about -40 kms -1 that is 
resolved only on some of the sight lines. The main component 
is probably associated with the Lupus molecular cloud region 
at 150 pc distance, while the second one may show gas in the 
Sco OBI shell. 

For the molecular gas we find strong arguments for a den- 
sity of the order of magnitude of 100 ctrT 3 . The radiation field 
Ua in the molecular cloud region may be strong, if one con- 
siders single-cloud H2 excitation models. The CO abundance 
rather hints to a radiation field near the mean Galactic value, 
which could also be consistent with multi-cloud models for the 
H2 excitation. At a high level of radiation in the clouds the 
molecular gas may be in a process of dissolution. 

In addition to this indicated dumpiness along the sight lines 
the gas structure perpendicular to the sight lines is analysed. 
The spatial variation in the column densities for the Lupus 
cloud gas is found to be 40% for H2 and 60% for Hi. This 
strongly indicates dumpiness in the gas on scales of 0. 1 pc and 
below. We are not able to resolve any spatial variation in the 
velocity of the absorbing gas over the field of view. 

dumpiness of the molecular interstellar medium has been 
detected in several recent works and on various scales, e.g., 
Pan et al. (2001) on sub-parsec scales and Lauroesch & Meyer 
(2000) on scales down to almost ten AU. 
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This work is another confirmation of the inhomogeneous 
structure of the Galactic ISM on sub-parsec scales. 

Appendix A: A simple model of H 2 excitation 

This Appendix gives a concise overview of the H2 excitation 
model used in our analysis (Sect. 4.2). The full details are found 
in Bluhm (2003). 

The basic properties of the model are: 

- 6 parameters (/V(H2), b, T, n p , n H , Ua) are used to calculate 
8 column densities N(J = 0..7). 

- Assumptions or approximations are: Ui(A) = const., 
7kin(0 = const., «h(0 = const., and n p (l) = const, (where I 
is the depth along the sight line in the cloud). 

- Collisional excitation by particles other than neutral and 
ionised hydrogen is neglected. 

- Collisional excitation and de-excitation of vibrationally ex- 
cited levels as well as UV absorption out of vibrationally 
excited levels are neglected. 

- The radiation which is responsible for the UV-pumping 
is assumed to come from the direction of the background 
source, i.e. the radiation propagates along the sight line. 

- While self-shielding in absorption lines is included (by in- 
tegration over the Voigt absorption profiles), blending of 
lines from different species (different /-levels of H2 and/or 
atomic species) is neglected. Thus for high column density 
sight lines, the radiation field will be somewhat underesti- 
mated. 

- Absorption of the UV continuum within the cloud is taken 
into consideration by applying approximate correction fac- 
tors to the cascade entry rates. The value of Ua refers to the 
irradiated edge of the cloud. 

- The UV-pumping process is calculated in a simplified man- 
ner by using averaged "cascade efficiency factors" (based 
on the cascade efficiency factors a(voJo; J) published in 
Black & Dalgarno 1976). The averaged factors a(J'; J) give 
the probability that a molecule in level J' before the ab- 
sorption of an UV photon ends up in level J after the UV- 
pumping process. 

- It is assumed that 87% of the absorptions lead to UV- 
pumping while the remaining 13% lead to dissociation (see 
Abgrall et al. 1992). In equilibrium, the total number of 
newly formed molecules is the same as the number of de- 
stroyed molecules, but the ortho-to-para ratio will differ in 
general. According to theoretical calculations, an ortho-to- 
para ratio of 3 is probable (see, e.g., Takahashi 2001). Since 
it is unlikely that a newly formed molecule has lost its for- 
mation energy completely when it leaves the dust surface, it 
will enter the rovibrational cascade, albeit with a different 
cascade entry distribution than the UV-excited molecules. 
Takahashi & Uehara (2001) modelled the vibrational and 
rotational excitation of H2 newly formed on icy, silicate, 
and carbonaceous dust surfaces. The most populated vibra- 
tional levels are v = 7, 6, and 2, respectively, with an aver- 
age rotational excitation of about J = 10. As cascade effi- 
ciency factors a j( J) for newly formed molecules mean val- 
ues of a(vo, Jo; J) (for vo = 2, 6, 7 and Jo = 9, 10) have been 



chosen. In an analysis of experimental results about the for- 
mation of molecular hydrogen Katz et al. (1999) find a non- 
vanishing probability that a newly formed molecule is not 
immediately ejected from the dust surface. A molecule lin- 
gering on the surface will loose a major fraction of its for- 
mation energy or even be "thermalized" before it is eventu- 
ally desorbed. To take this into account, for 40 % of the 
newly formed molecules a Boltzmann distribution at the 
gas temperature is assumed. 

Probably the most uncertain point of the model is the treatment 
of molecule formation. The rotational levels predominantly af- 
fected by this are J = 6 and 7 = 7. 
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